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ADAPTIVE  FILTERING  FOR  SEISMIC  SIGNAL  DETECTION 


Final  Report  (Part  I) 

I.  Statement  of  the  Problem 

This  project  is  cc.cerned  with  the  adaptive  digital  processing  to 
remove  (or  reduce)  noisy  signals  generated  by  sensors  at  large  distance 
away  from  the  sensors.  The  noisy  signals  need  to  be  removed  because  only 
those  signsals  generated  near  the  sensors  are  of  interest  for  detection. 

The  signals  to  be  detected  may  correspond  to  footsteps  of  a  walking  man.  - 
The  practical  approach  to  this  problem  at  the  Waterways  Experiment  Station 
is  to  use  two  or  more  spatially  separated  sensors.  Each  sensor  may  consist 
of  a  number  of  transducer  elements-each  contributing  constructively  to  the 
total  response  of  the  sensor.  Figure  A  is  a  sketch  of  the  sensor  layout. 

In  this  sketch,  the  seismic  energy  from  the  far-field  source  is  very  simi¬ 
lar  at  both  sensors  1  and  2.  Therefore,  we  should  be  able  to  reduce  the 
final  signal  used  for  processing  by  removing  the  correlated  portions  of 
the  signal  observed  at  sensors  1  and  2.  When  a  person  is  very  near  to 
sensor  1  (sensor  2  being  3  meters  away  or  more) ,  he  generates  a  much 
different  signal  in  sensor  1  than  in  sensor  2.  Therefore,  the  signals  of 
interest  should  not  be  correlated.  Mathematically  the  problem  is  to 
detect  the  footstep  signal  from  the  correlated  background  noise  caused 
by  e.g.  moving  vehicles  at  some  distance  away  plus  additive  sensor  noises. 
The  adaptive  procedure  is  needed  because  of  the  non-stationary  nature  of 
the  data  considered. 

II.  Methods  of  Approach 

Two  classes  of  digital  processing  techniques  are  considered;  namely 


the  adaptive  Kalman  filtering  and  the  adaptive  digital  filtering.  The 


adaptive  Kalman  filter  continuously  monitors  the  "state  jumps"  due  to 


the  presence  of  personnel.  Without  state  jumps  the  filter  performs  as 
an  ordinary  Kalman  filter.  The  filter  parameters  are  adjusted  if  a 
state  jump  is  detected.  False  detection  can  be  minimized  by  averaging 
over  a  sufficiently  large  "window.”  A  generalized  likelihood  ratio  is 
used  for  detection  of  state  jumps.  The  method  was  originally  proposed 
by  Willsky  and  Jones  [1]  and  was  successfully  employed  by  us  in  detecting 
the  object  boundary  in  image  scan  lines  [2]. 


The  block  diagram  of  the  adaptive  digital  filter  considered  is 
shown  in  Figure  B.  The  input-output  relation  of  the  non-recursive 
Adaptive  Digital  Filter  (ADP)  is  given  by  (Figure  C) , 

N 

g  -  l  b  f  (1) 

m  ~  n,m  m-n 

n=0 

where  g  is  the  filter  output,  b  is  the  filter  coefficient,  and  f  is 
m  n,m  n 

the  filter  input.  Define  the  vectors 


where  N  is  the  order  of  the  filter,  and  m  is  the  number  of  iterations  in 
adjusting  the  filter  coefficients.  Eq.  (1)  can  be  written  as 


g  -  B  '  F  =  F  '  B  (2) 

m  mm  mm 

We  use  the  Noisy  IMS  Algorithm  [3]  in  our  computer  program  to  calculate 


t. 
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where  v  is  a  constant.  After  adjusting  the  values  of  v,  N,  m,  and  delay 

time,  we  can  obtain  the  desired  result  with  a  greatly  improved  signal-to- 

noise  ratio. 

III.  Implementation  and  Computer  Results 

Extensive  effort  was  made  to  implement  the  two  methods  described 
above  on  the  PDP  11/45  computer  with  the  seismic  data  supplied  by  the 
Waterways  Experiment  Station.  The  data  are  stored  in  two  magnetic  tapes. 
The  first  tape  consists  of  4-channel  data  while  the  second  tape  consists 
of  14-channel  data.  Mr.  Mark  Flohr  of  the  Waterways  Experiment  Station 
was  particularly  helpful  to  acquaint  us  with  the  structure  of  the  data 
format  during  his  visit  in  August  1979.  The  computer  results  presented 
are  based  on  a  data  set  of  footstep  signal  plus  a  small  amount  of  noise. 

To  cope  with  the  memory  capacity  of  the  minicomputer ,  every  ’fifth  data 

sample  is  taken  from  the  4-channel  data  tape  with  a  resulting  sampling 
rate  of  200  samplej/sec.  A  sequence  of  960  sample  points  are  taken  for 
detailed  study  {Record  327  to  356).  Figures  la,  2a,  3a  and  4a  shew  the 
digitized  data  from  the  four  channels  with  corresponding  Fourier  ampli¬ 
tude  spectra  given  by  Figures  lb,  2b,  3b  and  4b,  respectively. 

The  Adaptive  Kalman  Filtering  was  implemented  to  detect  the  jumps. 

« 

Figures  5,6,  7  and  8  show  the  waveforms  before  ahd  after  detection  along 
with  a  printout  of  the  number  of  jumps  detected.  These  numbers  are  gene- 
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rally  correct.  The  complete  listing  of  the  computer  program  has  been 
forwarded  to  the  Waterways  Experiment  Station.  Further  improvement  of 
the  program  is  desirable  to  demonstrate  the  effect  of  noise  reduction 
while  performing  the  jump  detection.  The  parameters  are  fairly  sensi¬ 
tive  to  the  result  and  the  following  list  is  the  best  set  of  parameters: 
AA  (transition  coefficient)  *  1.0 
Q  (covariance  of  the  system  driving  noise)  =  1.0 
R  (noise  pcwer)  =  1.0 
TH  (threshold)  =  10.0 
M  (length  of  window)  =  10 

There  is  no  estimation  of  parameter  R  or  transition  matrices  W  and  P. 
Total  number  of  samples  is  960  per  channel. 

For  the  inqplementation  of  the  Adaptive  Digital  Filter  on  the  PDP 
11/45,  the  subroutine  AF  serves  as  the  Adaptive  Filter,  and  the  subrou¬ 
tine  SHIFTR  serves  as  the  Shift-Right  Register.  The  meanings  of  all 
variables  and  arrays  in  the  subroutine  AF  are  stated  below. 


ND: 

Noisy  data  array,  i.e.  input  data 

NF: 

Number  of  noisy  data  points 

NDELAY: 

Number  of  delay  points 

M: 

Number  of  iterations 

N: 

The  order  of  the  ADF 

V: 

The  constant  v  of  the  ADF 

B: 

Coefficients  of  the  ADF,  stored  in 

an  array 

DND: 

Delay  elements  for  storing  delayed 

data  points 

F: 

The  variables  f  in  the  ADF,  as  an 
m 

array 

G: 

The  variables  g  in  the  ADF,  i.e., 
m 

the  output 

E: 

The  variables  e  in  the  ADF 
m 

Y: 

The  variable  y 
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Array  B  and  array  F  are  initially  set  to  zero.  The  program  processes 
the  data  on  a  channel-by-channel  basis.  The  complete  program  named  as 
SPAN4  also  includes  an  ANFFT  Subroutine  which  performs  FFT  on  the  origi¬ 
nal  input  data  and  also  on  the  filtered  output  data.  The  complete  pro¬ 
gram  listing  is  in  the  Appendix. 

Figure  9  shows  a  sequence  of  filtered  waveforms  for  channel  1.  The 
best  result  is  given  by  Figure  9e  with  a  corresponding  FFT  given  by  Figure 
9f .  Here  M  =  1,  NDELAY  =  30,  N  =  18,  and  V  =  0.0112.  Figure  10  shows  a 
sequence  of  filtered  waveforms  for  channel  2.  The  best  result  is  given 
by  Figure  10e  with  a  corresponding  FFT  given  by  Figure  lOf.  Here  M  *=  1, 
NDELAY  =  30,  N  *  17  and  V  =  0.00145.  It  is  noted  that  the  parameter 
choice  is  also  very  sensitive  in  the  adaptive  digital  filtering  operations. 
The  filtered  results  clearly  illustrate  the  suppression  of  background  noise 
and  the  enhancement  of  footstep  signals.  The  filtered  results  are  so  good 
that  further  process  of  jump  detection  appears  unnecessary.  All  Fourier 
spectra  correspond  to  a  frequency  range  of  0  to  100  Hz. 

Both  adaptive  filtering  techniques  have  been  shown  to  be  very  effective 
for  the  seismic  detection  tasks  under  consideration.  They  both  are  computa¬ 
tionally  efficient.  Proper  choice  of  parameters  is  an  important  step  to 
derive  the  best  results.  Although  mathematically  the  Kalman  filter  and 
the  adaptive  digital  filter  are  somewhat  equivalent,  the  implementations 
of  the  two  techniques  are  different  and  thus  the  results  are  different  also. 
The  adaptive  Kalman  filtering  in  fact  has  more  flexibility  to  adjust  itself 
to  the  change  in  data  characteristics. 
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END 
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